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Abstract 

Within the framework of the standard model we observe that there is 
a significant discrepancy between the most precise Z boson decay asym- 
metry measurement and the limit from direct searches for Higgs boson 
production. Using methods inspired by the Particle Data Group we ex- 
plore the possible effect on fits of the Higgs boson mass. In each case 
the central value and the 95% confidence level upper limit increase sig- 
nificantly relative to the conventional fit. The results suggest caution in 
drawing conclusions about the Higgs boson mass from the existing data. 
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Introduction Perhaps the most pressing question in particle physics today 
is the mass scale of the quanta that break electroweak symmetry, giving mass 
to the particles in the theory, including the quark and lepton constituents of 
ordinary atomic matter. That scale determines whether the symmetry breaking 
force is weak or strong and it sets the energy scale future accelerators will need 
for detailed studies of the mass-generating mechanism. In general the issue can 
only be resolved by discovering the symmetry breaking quanta at a high energy 
collider. However, in particular theoretical frameworks, such as for instance 
the standard model, radiative corrections to already measured quantities can be 
used to constrain the mass of the symmetry breaking sector. 

Interpreted in the standard model framework, beautiful data from LEP, 
SLAC, and Fermilab appear to favor a light Higgs boson with mass of order 
100 GeV.0 The conclusion emerges from the effect of virtual Higgs bosons, via 
radiative corrections, on precision measurements of the Z and W bosons. In 
addition, the four LEP experiments have searched for real Higgs bosons, with 
negative results that when combined are expected to imply a lower limit run > 
77 GeV at 95% confidence level. @ Taken together the experiments suggest a 
window between 80 and a few hundred GeV. The purpose of this note is to 
suggest that the window may in fact be substantially larger, in part because of 
well known inconsistencies within the precision data, but more because of equally 
significant inconsistencies between precision data and the direct searches whose 
magnitude has, with some noteworthy exceptions, || || gone largely unnoticed 
and/or unremarked. 

The problem of how to combine inconsistent data has led to the break-up 
of many beautiful friendships. The mathematical theory of statistics provides 
no magic bullets and ultimately the discrepancies can only be resolved by future 
experiments. The Particle Data Group f| (PDG) has for many years scaled 
the uncertainty of discrepant results by a factor I will call Spdg 5 defined by 
Spdg = \Jx 2 1 [N — 1), where N is the number of data points being combined. 
They scale the uncertainty of the combined fit by the factor Spdg if and only 
if Spdg > 1- This is a conservative prescription, which amounts to requiring 
that the fit have a good confidence level, ranging from 32% for N = 2 to greater 
than 40% for larger values of N. If the confidence level is already good, the 
scale factor has little effect; it only has a major effect on very discrepant data. 
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The PDG argues (see |J) that low confidence level fits occur historically at 
a rate significantly greater than expected by chance, that major discrepancies 
are often, with time, found to result from underestimated systematic effects, 
and that the scaled error provides a more cautious interpretation of the data. 
A few authors |3|, [7|, |^| have applied Spdg to the asymmetry measurements of 
sin 2 8 l eff° n , as I will also do here. 

With the top quark mass fixed at the value determined by CDF and DO, the 
most sensitive probe of m# is currently the effective leptonic weak interaction 
mixing angle, sin 2 6 l ^Jj° n , measured in a variety of Z boson decay asymmetries. 
The extent to which the asymmetries currently dominate the estimate of m# can 
be seen by comparing the conventional (unsealed) fit to the seven asymmetry 
measurements, mu = 104±54°GeV, with the LEP electroweak working group 
global fit to all data, m H = 115±gg 6 GeV. Because it more than suffices for the 
purposes of this paper, the analysis that follows is based on the asymmetry 
measurements alone. The PDG scale factor then increases the uncertainty but 
not the central value of the combined fit for sin 2 9 l ^ on and m#.f] 

The focus of this paper is on the discrepancies between precision measure- 
ments and the limit from the direct searches, which will be addressed by a 
method analogous to the PDG scale factor. Like the PDG prescription, the 
idea is to scale the error so that the precision measurement has a significant 
probability P to be consistent with the direct search limit. I consider P = 0.32, 
corresponding to the PDG's choice, as well as larger and smaller values. To ac- 
count for uncertainty in the search limits, which may also be subject to unknown 
systematic errors, I consider a range of different lower limits on mn, from a very 
conservative 50 GeV to a futuristic 90 GeV. In this approach both the central 
value and the uncertainty of the fit are affected. In addition I present fits using 
two other methods discussed by the PDG. 

By giving full weight to a measurement that is in serious conflict with the 
direct search lower limit, the conventional method risks underestimating m#. 
The alternative methods considered here provide a more conservative estimate 
of the upper limit on m# but risk skewing the fit to large m#. Taken together 

3 If instead the rescaled fit to sin 2 8 l ^ on is included in a global fit to all Z decay data, as 
in Gurtu's analysis the rescaling does affect the central value of mn . 
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the results strongly suggest caution in drawing conclusions from the precision 
data about the value of the Higgs boson mass. 

The precision data The relevant values of sin 2 9 l ^ on and the quoted experi- 
mental uncertainties are shown in table 1, from the preliminary values presented 
at the 1997 summer conferences. |l|] For each value the table displays the corre- 
sponding value of ran and the 95% CL upper (m^) and lower (m^ 5 ) bounds (that 
is, the symmetric 90% confidence intervals). Also indicated is the probability 
for rriH to lie below 77 GeV. Gaussian distributions are assumed for sin 2 6 l ^ on 
and log{mn)- 



The values of m# are from the state of the art MS computation of reference 



H - see also |T(J. To obtain the confidence intervals and probabilities the 
parametric error is combined in quadrature with the experimental errors. The 
parametric error is equivalent to ±0.00030 uncertainty in sin 2 6 l ^ on — - see 
|8|, 10 1 . It is dominated by roughly equal contributions from the uncertainties 



in the top quark mass, m t = 175 ± 6 GeV, and the fine structure constant 
at the Z mass, a~ l (mz) = 128.896 ± 0.090, in addition to other much smaller 
contributions, including A«qcd(^z) and uncomputed higher order corrections.^ 
The six LEP measurements in table 1 are each combined from the four LEP 
experiments, and in each case the combined fit has a good confidence level. The 
conventional maximum likelihood fit for the LEP measurements is shown in the 
first row of table 2. The chi-squared per degree of freedom is y 2 /N — 1 = 4.4/5 
corresponding to a robust 0.5 confidence level. The central value is mu = 240 
GeV and the 95% CL upper limit is 860 GeV. There is no entry for the Spdg 
fit since Spdg < 1- 

Combining all seven measurements (the conventional LEP + SLC fit in table 
2) the central value decreases to 100 GeV and the 95% CL upper limit falls to 310 
GeV, demonstrating the effect of the high precision and lower sir?Q l ^ on from 
Alr- The chi-squared per degree of freedom now rises to 12.5/6, with a marginal 



4 There are also negligible extrapolation errors from reference B, equivalent to < 0.00003 
in sin 2 l ^ff On for 75 < tuh < 600 GeV. Even outside this range they have no real effect on the 
analysis, since the confidence levels and scale factors only depend on the relationship between 
sin 2 9 l ^ on and mu for mu = m 1 ^ 111 *. The worst case is then m x ™ lt = 50 GeV, close enough 
for any additional error to be negligible. The very large values of mg" 5 in the tables could be 
affected but they have no precise significance in any case. 
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confidence level of 0.05. The PDG scale factor is then Spdg — 1-45. Using it, 
the combined uncertainty of the fit increases from ±0.00023 to ±0.00033 and 
the 95% CL upper limit on m H increases modestly to 420 GeV. 

. . . and the direct search limits In addition to discrepancies among the mea- 
surements of sin 2 9 l eff° n , which are problematic whether we assume the standard 
model or not, table 1 also reveals a second discrepancy that occurs specifically 
within the standard model framework. The most precise measurement, Am, 
implies a 95% CL upper limit on m# of 77 GeV, while the direct searches at 
LEP are expected to provide a combined 95% lower limit also at 77 GeV. 
(The individual 95% CL limits quoted by the four experiments range from 66 
to 71 GeV.0) The third most precise measurement, A l FB , also has significant 
weight (71%) below the direct search limit.f] 

This raises a difficult question: within the standard model framework what 
role if any should the direct search limits play in extracting the implications 
of the precision data? There is no single "right" answer to the question. A 
maximum likelihood fit including both the precision data and the direct search 
data would replicate the conventional fit if the central value lies above the lower 
limit, m 1 ^™ 1 *, from the direct searches. That is a defensible interpretation, since 
if the true value of m# were near m l ™ lt we would expect values of run obtained 
from measurements of sin 2 6 l ^Jj° n to lie both above and below m 1 ^ 111 *. By under- 
weighting downward fluctuations while leaving upward fluctuations at their full 
weight, we risk skewing the fit upward. Mindful of this risk, it is still instructive 
to explore the sensitivity of the fit to the weight ascribed to measurements that 
are individually in significant contradiction with the direct search limit. 

Clearly the direct search results are not irrelevant. If, for instance, the only 
information available were the direct search limits and the Am measurement, 
we would conclude that the standard model is excluded at 90% CL.[] Theorists 
would have flooded the Los Alamos server with papers on the death of the 
standard model and the birth of new theories W,X,Y,Z... In the actual situation 
the Alr measurement causes the fit to run to shift by more than a factor two, 
from 240 GeV to 100 GeV, and the 95% upper limit to fall from the TeV scale 

5 Aii{ and A l FB are also the only measurements with mg" 5 below the TeV scale. 
6 I thank Lawrence Hall for this perspective. 
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to ~ 340 GeV. It is fully weighted in the conventional standard model fit despite 
a significant contradiction with the standard model. 

If the discrepancy were even greater — say, for instance, a precision mea- 
surement implying m# = 10 MeV with a 99.99% CL upper limit at 77 GeV - 
the clear response would be to omit that measurement from a standard model 
fit, although it could still be considered in a broader framework encompassing 
the possibility of new physics.[] On the other hand, A l FB , with 31% probability 
to be consistent with a 95% lower limit at 77 GeV, would surely be retained. 
The question is how to resolve the intermediate cases in which the discrepancy 
is significant but not so significant that the data should clearly be excluded. 

Consider a prescription, analogous to the PDG scale factor, that interpo- 
lates smoothly between the extremes. Imagine a measurement x with experi- 
mental error 5e and a quantity y that is related to x with an uncertainty x ± Sp 
(the parametric error). Suppose there is an upper limit on y at y = yo that trans- 
lates to an upper limit on x at x ± Sp, such that the measurement x falls below 
the implied limit, x < Xq. The discrepancy between the measurement and the 
limit is then characterized by a Gaussian distribution centered at x with stan- 
dard deviation a = \J S 2 E + Sp, with a computable probability P for x > xq. If 
P is less than a chosen minimal confidence level Pyp (VR for "virtual- real" ) , 
then 5p is scaled by a factor Svr chosen so that the Gaussian centered at x 
with standard deviation a' = J (SvrSe) 2 + Sp has probability Pvr for x > Xq. 
If xq — x is small enough, the scale factor has little or no effect. If x is many 
a below x , Svr will be large and the data point x will have reduced weight 
in a combined fit with other data. Intermediate cases will interpolate smoothly 
between the two extremes, depending on the values of x — xo, a, and Pvr- 

The value of Pvr is of course arbitrary. One plausible choice is Pvr = 0.32 
(or 0.3173 to the cognoscenti), since that is the confidence level implicit in the 
PDG scale factor for N = 2. A plausible choice for the lower limit on irtn is 
m hmit _ 77 GeV. The resulting fits are shown in table 2. The fit to the LEP 
data is affected only modestly, with an increase of ~ 10% in m#. For the LEP 
+ SLC fit, the central value of win and the 90% confidence interval increase 

7 In fact, parity violation in atomic Cesium currently implies m# ^11 MeV (MeV is not a 
typographical error) though only 1.2a from 77 GeV.|jll| Its weight in the combined fit would 
be negligible. 
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significantly, to the values of the LEP fit. The scale factors are Svr(A 1 fb ) = 1.2 
and S VR (A LR ) = 4.1. 

Table 3 displays the results of varying Pvr = 0.20, 0.32, 0.40 and m 1 ™ 1 * = 
50, 60, 70, 80, 90 GeV. Though not reported by the experimental groups, the 
confidence levels at 50 and 60 GeV are probably much tighter than 97%, both 
because the LEP II cross sections increase for smaller m# and because the LEP 
I data contributes to the confidence level at those masses — - the 95% CL lower 
bound from LEP I alone is > 66 GeV. The value 80 GeV is close to the 
presently projected 95% combined limit of the four LEP experiments, while 90 
GeV is the anticipated limit if no discovery emerges from currently planned 
LEP II running. Vacant entries in table 3 indicate that the fits are unmodified, 
Svr < 1) and that the conventional fit, line 1 of table 2, applies. For the 
LEP data, the Higgs mass scale varies by no more than a factor 1.5 from the 
conventional fit over the entire range of table 3. For the LEP + SLC data, the 
difference is a factor 1.5 for (Pvr, m}™ 1 ) = (0.20, 50 GeV) and becomes as large 

factor 4. 

Other methods In this section I will briefly present results using methods 
discussed by the PDG[|l2j for combining measurements that conflict with a limit. 
They are no less arbitrary than the Svr scale factor method discussed above. 

Consider a collection of measurements Xi ± 5e,i, i = l,...,iV, some of 
which are nominally inconsistent with an upper limit at Xq. The "Bayesian" 
method is to combine all data points in the conventional way and to multiply 
the combined Gaussian distribution by a step function C6(x — x ), so that the 
distribution vanishes below xq. C is a normalization factor to guarantee total 
unit probability. I have modified the usual Bayesian method to account for the 
fact that the lower limit is not absolute but has 95% confidence, by choosing C 
to give the distribution probability 0.95 for x > xq- The 50%- and 95%-tiles of 
the resulting distributions are shown in table 2 for m 1 ^ 1 * = 77 GeV. 

Three "frequentist" prescriptions are also discussed by the PDG, of which 
two are considered here. For x^ < x one prescription assigns Xjit^j — > x ±Se,i 
when the limit Xq is known exactly. Including the parametric uncertainty, I 
modify this to X{ ± 8e,{ — * xq ± &i where Oi = ^JS E i + 5 p. The readjusted 
points are then combined as usual (including Spdg if applicable) with the other 
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measurements. An extremely conservative variation, intended only to obtain 
the 95% CL upper limit, replaces Xi — > min(xj, x + 1.64<7j), so that 95% of the 
probability distribution for each measurement is above the limit x . The results 
are illustrated in table 2 for m^ mit = 77 GeV. 

Conclusion Several related points are deferred to a more detailed report, 
including the following: (1) There are indications, depending on how the data 
is grouped, that the confidence level of the conventional fit may be even less 
than 0.05. Dependence on how the data is grouped reflects the uncertainty 
in the confidence level for fits of small data samples. (2) Though definitive 
conclusions can only come from the experimental groups, two estimates suggest 
that improved 6-tagging methods are not likely to cause a big shift in A b FB . (3) 
The W boson mass measurement currently has a sensitivity to ran at roughly 
the middle of the pack of the asymmetry measurements. It does not qualitatively 
alter the conclusions. 

In summary, the A LR measurement is inconsistent at 95% CL both with 
the LEP asymmetry measurements and, in the standard model, with the Higgs 
boson search limits, while its precision causes it to have a profound effect on the 
combined standard model fit. The conflict with the search limits may diminish 
or disappear if there is new physics outside the standard model framework but 
is necessarily germane to a standard model fit. The analysis presented here 
is meant as a warning signal, a yellow if not a red flag, suggesting caution in 
drawing conclusions from the precision data about the mass of the standard 
model Higgs boson. Applying methods inspired by the Particle Data Group to 
these discrepancies, we find that the central value of m# increases by factors 
from ~ 1.5 to ~ 3 while the 95% CL upper limit increases toward the TeV scale. 
Only future experimental results can resolve the discrepancies in the present 
experimental situation. 
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Tables 



Table 1. Values for sin 2 9^j on from asymmetry measurements |1|] with la 
experimental errors. The corresponding Higgs boson masses, the 95% CL upper 
and lower limits, and the confidence level for m# < 77 GeV are given for each 
measurement. 





sin^r™ (1*) 


m H (GeV) 


m95> m 95 


P(< 77GeV) 


Alr 


0.23055 (41) 


16 


3, 80 


0.95 


A b 

FB 


0.23236 (43) 


520 


100, 2700 


0.03 


A 1 

FB 


0.23102 (56) 


40 


5, 290 


0.71 


A T 


0.23228 (81) 


440 


30, 6700 


0.14 


A e 


0.23243 (93) 


590 


28, 13000 


0.14 


Qfb 


0.23220 (100) 


380 


14, 10000 


0.21 


A c 

FB 


0.23140 (111) 


83 


2, 3000 


0.49 
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Table 2. Fits to the LEP and LEP + SLC data as described in the text. 
The Bayesian and frequentist fits assume m 1 ^ 111 * = 77 GeV. 



Data set 


Fit 


sin 2 ^ e ff pton (l(T) 


m H (GeV) 


m^g, m> 5 




Conventional 


0.23196 (28) 


240 


67, 860 




Spdg 








LEP 


S VR (0.32, 77) GeV) 


0.23203(29) 


270 


74,1000 




Bayes 


0.23197 


250 


< 880 (95%) 




Frequentist (1) 


0.23204(28) 


280 


78,1000 




rrequentist (2) 


< 0.23304 (95%) 




1900 (95%) 




Conventional 


0.23152 (23) 


100 


32, 340 




Spdg 


0.23152 (33) 


100 


26, 420 


LEP + SLC 


S VR (0.S2, 77) GeV) 


0.23198 (28) 


250 


69,920 




Bayes 


0.23171 


150 


< 500 (95%) 




Frequentist (1) 


0.23183 (23) 


190 


58,610 




Frequentist (2) 


< 0.23284 (95%) 




< 1300 (95%) 
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Table 3. Fits to LEP and LEP + SLC data using the S VR scale factor 
with various values of P V r and m 1 ™ 1 *. Each entry displays the central value of 
m H and m^, m^, the 95% CL lower and upper limits, in GeV. Empty entries 
indicate that no measurement is far enough below threshold to be modified by 
the scale factor and that the conventional fit of table 2 applies. 





LEP 


LEP + SLC 


™ limit 

m H 


P VR = 0.20 0.32 0.40 


0.20 0.32 0.40 


50 




160 210 230 
41,590 59,740 65, 830 


60 


310 
82, 1100 


170 220 300 
47, 650 61,770 80, 1100 


70 


360 
94, 1400 


190 220 350 
51, 680 62,790 91, 1300 


80 


290 380 
77,1100 100, 1500 


190 260 370 
54,700 72,970 96,1400 


90 


320 390 
85,1200 100, 1500 


200 300 380 
56, 710 79,1100 100, 1500 
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